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$—( \ Abstract. In two recent reactions by Belle producing DD and DD* meson pairs, peaks above 

MH' threshold have been measured in the differential cross sections, possibly indicating new resonances 

■^ ' in these channels. We want to study such reactions from the point of view that the D meson pairs 

are produced from already known or predicted resonances below threshold. Our study shows that 

the peak in the DD* production is not likely to be caused by the X(3872) resonance, but the peak 

\ seen in DD invariant mass can be well described if the DD pair comes from the already predicted 

I— I . scalar X (3700) resonance. 

(~) I Keywords: Dynamically generated resonances, X (3872) 

^ I. INTRODUCTION 

,__i \ Belle collaboration has measured reactions producing D mesons together with a 7/i// 

> ' from e+e collisions [1]. The invariant mass distributions for DD and DD* in these 

^ . reactions peak around 70 MeV above threshold. This observation lead to the claim 

^^ \ of new resonances. In a previous paper [2] Belle has studied the DD* invariant mass 

CN \ distribution in B decays {B -^ KDD*), in this case a peak was observed only a few 

Tj- ■ MeV above threshold and it was first suggested that this peak could be indicative of the 

O I existence of a new X (3875) resonance. 

Q I This second experiment has been analyzed in several theoretical works. In [3] it is 

l" • concluded that the peak seen in [2] may be caused by the X (3872) if it is a virtual state. 

• ^ . A latter work [4] shows that this conclusion may be changed if one takes into account 

rS I explicitly the D* width. A different approach is suggested in [5], where the X(3872) is 

c^ ■ generated dynamically by the interaction of pseudoscalar mesons with vector ones. 

The idea of generating dynamically resonances in order to describe the newly found 
charmed states is not new and had already been exploited in [6, 7, 8, 9]. We have 
extended the ideas of these works including the hidden charm sector in [5, 10]. 

This work presents the results of [1 1] in which the peaks seen in [1] are interpreted as 
caused by some of the hidden charm subthreshold resonances dynamically generated in 
[5, 10]. 

II. FRAMEWORK 

We assume that the D mesons in the reaction e+e -^ J /xifDD come from a resonance. 
In this case the reaction can be described by the diagram in figure 1 and close to the DD 



threshold the amplitude depends strongly only on the DD invariant mass. 
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FIGURE 1. Feynman diagram for the process e+e^ -^ J/xj/DD 

This amplitude can be written in terms of the X propagator with all other parts 
factorized as a constant C: 



^ = c = CD 

MfjDD)-Ml + irxMx' 

if we describe the X resonance as a Breit-Wigner type. 

The cross section can be calculated from the amplitude by integrating over the phase 
space of the three particles in the final state. Since ^ depends only on the DD invariant 
mass, one can evaluate the differential cross section that takes the form: 

'" - ' "-l^ll7?ll^l^ (2) 



dMi„^{DD) {2nys^' 

where s is the center of mass energy of the electron positron pair squared and | k \ and 
1 7? I are given by the relative momentum of the D mesons and the relative momentum of 
the 7/ 1// and the D meson system in the rest frame of each system respectively. 

Now we are going to describe the model for generating dynamically the resonances. 
First a field is written for all mesons belonging to a 15-plet of SU{4), one field, $, for 
the pseudoscalars and another, ^ for the vector mesons. Explicit expressions for these 
fields can be found in [5, 10]. 

The Lagrangians are built by coupling two currents of the meson fields. The La- 
grangians built this way are SU{4) invariant, but since this is not an exact symmetry 
in Nature, one breaks it by suppressing the terms of these Lagrangians where the under- 
lying interaction is driven by the exchange of a heavy vector mesons [12, 13]. 

One collects the tree level amplitudes for meson-meson scattering, projected in s- 
wave, for each couple channel space, into a matrix. This matrix is used as input to 
solve the Bathe- Salpeter equation that in the on-shell formalism of [15, 16], assumes 
an algebraic form: 

T = (i -yG)"V, for scalars (3) 

T = -(i+yG)"V"e'."e'',foraxials (4) 

The loop function, Gu, has its imaginary part fixed in order to ensure unitarity of 
the T-matrix. Resonances are identified as poles in the second Riemann sheet of the 
T-matrix. 



TABLE 1. Results of Mx and x for different values of an- 
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Within this model the D*q(23 17) , D^i (2460) and X(3872) are well reproduced. More- 
over interesting predictions are a heavy hidden charm scalar with mass around the DD 
threshold and an axial state with negative C-parity, almost degenerated in mass with the 
X(3872). Also broad open charm states with exotic quantum numbers appear. These re- 
sults are, though, less stable and contrast with previous works done where states with 
the same quantum numbers are narrower [6, 7, 8, 9]. 

The way we are going to apply this model to describe the reaction of figure 1 is to 
substitute the Breit-Wigner amplitude of eq. (1) by the DD T-matrix calculated from eq. 
(3). For the reaction e+e^ -^ J/\ifDD* everything is done analogously using eq. (4) for 
the T-matrix. 



III. RESULTS 



The e+e process measured by Belle [1] has a center of mass energy of \/s — 10.6 GeV. 
The Minv{DD^*') is measured from threshold up to 5.0 GeV. Since our model is reliable 
for energies within a few hundreds of MeV around the thresholds, we compare our 
model with belle's data for the points below 4.2 GeV. The invariant mass distributions 
are measured in counts per bin. To simulate the bins we integrate our theoretical curves 
in bins of the same size as experiment and normalize so that the total integral of our 
curve matches the total number of events in the experiment. 

Finally we compare the data with our predictions by a standard X^ test. 

The model has one free parameter, which is the subtraction constant in the loop 
function, «//. We show results for different values of an which correspond to a natural 
size. For the case of two D mesons in a loop, for example, the value aH=- 1 .3 correspond 
to a cut off with q^ax = 850 MeV. We present results for different values of an together 
with the corresponding pole position. 

In table 1 we show results, in the scalar sector, for the value of x^ calculated with the 
data from Belle, for all points below 4.2 GeV in the J /x^rDD and J/\irDD* production. 
Fig. 2 shows plots of our theoretical histograms compared with the experimental data 

[1]. 

One can see that the values of x^ arc around 1 for the X(3700), indicating a good 
fit to the data. Although the large experimental errors are in part responsible for this 
good fit the data is compatible with the existence of the scalar resonance below the DD 
threshold. 

In the case with a DD* pair in the final state we have to apply our model for the 
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FIGURE 2. Left: Theoretical histograms compared with data from [1] for DD invariant mass distribu- 
tion. Right: Theoretical histograms compared with data from [1] for DD* invariant mass distribution. 



interaction of pseudoscalars with vector mesons, in this case the resonances generated 
are axials, and the X in figure 1 should be identified with the X(3872) generated 
dynamically in our model. Table 1 also shows results for Mx and x^ for different values 
of Uh for this sector. Since the state X(3872) is known and has a rather precise mass, 
we have chosen a smaller range to vary the parameter a. 

In this case the x^ obtained is in all cases bigger than 3, clearly indicating a poor fit 
to the data. 

In the experimental paper the peaks have been fitted with Breit-Wigner type reso- 
nances [1], suggesting two new States. We have also performed such a fit in order to com- 
pare the results and have more meaningful conclusions. We take the same Breit-Wigner 
parameters suggested in the experimental paper. The scalar resonance with Mx=3878 
MeV and r;^=347 MeV and the axial one with Mx= 3942 MeV and Yx= 37 MeV. We 
show the results obtained by fitting a Breit-Wigner form from eq. (1) in ^ of eq. (2) in 
figure 3. Additionally we calculate x^ ^^d find x^/d.o.f =2.10 for the DD distribution 
and x^/d-o.f =1.34 for the DD* distribution. The value of X^ for the DD distribution 
can be improved if we take different parameters for the Breit-Wigner resonance. Taking 
for the fit Mx=3750 MeV and rx=250 MeV we obtain a value of x^/d.o.f =1.12, still 
slightly bigger than those obtained in our previous analysis. The value of x^ for the DD* 
distribution is undoubtedly better in the case of a Breit-Wigner fit that in our analysis 
assuming the X (3872) resonance as the X in the mechanism of fig. 1. 
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FIGURE 3. Left: Histograms calculated with Breit-Wigner resonance with mass Mx=3940 MeV com- 
pared to data. Right: Histograms calculated with Breit-Wigner resonance with mass Mx=3880 MeV 
compared to data. 



IV. SUMMARY 

We have analyzed a recent experiment from Belle [1] where peaks are observed in the D 
meson pair invariant mass suggesting new states with hidden charm. We have calculated 
the differential cross sections for these measurements from the point of view that the D 
meson pairs are produced by dynamically generated resonances below threshold. 

Our theoretical model describes very well the data for the DD invariant mass distribu- 
tion as being caused by a dynamically generated scalar resonance below DD threshold. 
On the other hand, the X^ values obtained in the fit of our theoretical curves for the DD* 
invariant mass distribution are too big, favoring thus the interpretation of a new state 
with mass around 3940 MeV. 

Although we have a good description of the scalar resonance, the large experimental 
uncertainties do not discard other possibilities. New experiments to search for this 
predicted scalar state are most welcome. 
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